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I. INTRODUCTION 



Matter enhanced oscillations of neutrinos via the 
Mikheyev, Smirnov, Wolfenstein (MSW) mechanism 0] 
is of considerable current interest. MSW mechanism may 
account for the deficit of solar neutrinos Q. Matter- 
enhanced neutrino oscillations may also play an im- 
portant role in neutrino propagation through the core- 
collapse supernovae 

The equations of motion for the neutrinos in the MSW 
problem can be solved by direct numerical integration, 
which must be repeated many times when a broad range 
of mixing parameters are considered. This often is not 
very convenient; consequently various approximations 
are widely used. Exact or approximate analytic results 
allow a greater understanding of the effects of param- 
eter changes. Exact solutions of the MSW equations 
for a number of density profiles are available in the lit- 
erature. Landau-Zener method [Q is applicable to the 
level-crossing phenomenon inherent in matter-enhanced 
neutrino oscillations (^J^. The standard Landau-Zener 
method is exact a linear density profile. Exact solutions 
for an exponential density profile or for a profile in 
the form tanha; pO| are also given in the literature. 

The analogy between supersymmetric quantum me- 
chanics and matter-enhanced neutrino oscillations was 



pointed out in Ref. |Tl|. Indeed this analogy was used 
to obtain a semiclassical approximate expression for the 
hopping probability ||ll|,|l2| . The aim of the current paper 
is to exploit this analogy further to obtain exact solutions 
for a class of electron density profiles. The integrability 
condition we investigate is analogous to the integrability 
condition called shape-invariance in supersymmetric 
quantum mechanics. 

for the purpose of establishing notation the MSW effect 
is outlined and salient formulas are given in Section II. 
In Section III the consequences of the shape-invariance of 
electron density profile are discussed, an exact expression 
for the electron neutrino amplitude for such density pro- 
files is derived and its asymptotic limits are calculated. 
Several examples are given in Section IV and existing 
results in the literature are shown to follow from the for- 
malism developed in the previous section. In Section V, 
the shape-invariance condition is shown to be a suitable 
first step for obtaining approximate expressions for the 
neutrino survival probability. Section VI includes a brief 
discussion of the results. 



II. OUTLINE OF THE MSW EFFECT 

The evolution of flavor eigenstates in matter is gov- 
erned by the equation 
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where we defined 
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In these equations Srr? = — ra\ is the vacuum mass- 
squared splitting, 9v is the vacuum mixing angle, Gf is 
the Fermi constant and Ne{x) is the number density of 
electrons in the medium. 

By making a change of basis 



, (2.4) 
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sin6'(a;) cos6'(a;) 







the flavor-basis Hamiltonian of Eq. (2.1) can be instanta- 
neously diagonalized. The matter mixing angle is defined 
via 



and 
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In the adiabatic basis the evolution equation takes the 
form 
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(2.7) 



where prime denotes derivative with respect to x. Since 
the 2x2 "Hamiltonian" in Eq. (2/7) is an element of the 
SU{2) algebra, the resulting time-evolution operator is 
an element of the SU{2) group. Hence it can be written 
in the form 



U = 



^i{x) -^*^{x) 

^2{x) nix) 



(2.8) 



where 5'i(x) and ^'2(2;) are solutions of Eq. (2.7) with 
the initial conditions ^'1(0;^) = 1 and '^2{xi) = 0. 

To calc ulate the electron neutrino survival probability 
Eq. (^.iD needs to be solved with the initial conditions 



1 and = 0. Using Eq. (2.8) the general solution 



(2.1) satisfying these initial conditions can be written as 



^'e(x) = cos6'(a;)[cos6'i^'i(a;) - sin6'i^'2(x)] 

-I- sin6l(x)[cos6',*2(a;) +sin6',^'t(a;)], (2.9) 

where 9i is the initial matter angle. Once the neutrinos 
leave the dense matter (e.g. the Sun), the solutions of 
Eq. (2.7) are particularly simple. Inserting these into 
Eq. ( |2.9D we obtain the electron neutrino amplitude at a 
distance L from the solar surface to be 



^'e(L) = cos 
+ sin6',„[cos 6^1*2, (S) 



COs6',*i^(5.) 

sin 6'j*i_(5)] exp 



exp I i , _ L 
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(2.10) 



where ^'i^(S) and ^2,(S) are the values oi'^i{x) and ^'2(x) 
on the solar surface. The electron neutrino survival prob- 
ability averaged over the detector position, i, is then 
given by 



P {u,^v,) = {\^,{L)\^)l 
1 1 

= 2 + 2 
1 



cos2e',„cos26', (l-2|«'2,(s)n 



■ cos 29„ sin 26',- 



(*l,(S)*2,(S) + *t^(S)*2,(S)) • (2-11) 



If the initial density is rather large, then cos 26 i ^ — 1 
and sm.29i ^ and the last term in Eq. (2.11) is very 
small. Different neutrinos arriving the detector carry 
different phases if they are produced over an extended 
source. Even if the initial matter density is not very 
large, averaging over the source position makes the last 
term very small as these phases average to zero. The 
completely averaged result for the electron neutrino sur- 
vival probability is then given by 



P{ue ^i^) ^ ^ + ^ cos29y{cos29, 



where the hopping probability is 

Phop = |*2,(S)P, 



I /source (1 2_P[^op) 5 

(2.12) 
(2.13) 



obtained by solving Eq. ( p.7[ ) with the initial conditions 
"ifiixi) = 1 and 1'2(a;i) = 0. 



III. SHAPE-INVARIANCE OF THE ELECTRON 
DENSITY PROFILE 

The coupled first-order equations for the flavor-basis 
wave functions can be decoupled to yield a second order 
equation for only the electron neutrino propagation 
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-[A + ip^{x) +ihip'{x)\'^e{x) =Q. (3.1) 
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Introducing the operators 
d_ 

dx 
d 



yl_ = ih- 



ih- 



dx 



Eq. {2_1) takes the form 
i-*e(a;) = 



fix) , 



A+-^f,{x) = %/A^'e(a;), 



and Eq. (3.1) can be compactly expressed as 



(3.2) 



(3.3) 



(3.4) 



In these equations we take (p(x),A^, and A- a function 
of a set of parameters a that specify space-independent 
properties of the electron density. 

The analogy between Eq. ( |3.lD and supersymmetric 
quantum mechanics was pointed out some time ago . 
In a paraUel development, in the context of supersym- 
metric quantum mechanics it was shown that a subset 
of the potentials for which the Schrodinger equations are 
exactly solvable share an integrability condition called 
shape invariance 0. In this paper we consider a gen- 
eralization of the shape-invariance condition to electron 
densities. We call a given electron density to be shape- 
invariant if a change of parameters satisfies the condition 

— f'^{x; ai) + ih(p'{x; ai) 

= -ip^{x] 02) - ihip'{x] 02) + R{ai) , (3.5) 

where 02 is a function of oi, and the remainder R{ai) 
is independent of x. As we illustrate in the forthcom- 
ing sections a number of electron density profiles satisfy 
this condition. Here we take oi to be the original set 
of parameters specifying the given electron density pro- 
file (hence real). The new set of parameters 02 are, in 
gene ral, complex. The shape invariance condition of Eq. 
(3.5) ca n b e rewritten in terms of the operators defined 
inEq. (pi) 



i_(ai)i+(ai) = A+{a2)A^{a2)+R{ai). 



(3.6) 



We assume that replacing oi with 02 in a given operator 
can be achieved with a similarity transformation: 

r(ai)0(ai)f-i(ai) = 0(a2). (3.7) 

If the parameters ai and a2 are related by a translation, 

(3.8) 



02 = ai 



then the the operator T(ai) of Eq. (3.7) is simply given 
by 



T{ai) = exp 77 



d_ 



(3.9) 



One should emphasize that the discussion in this article 
is valid for any relationship betwe en a i and 02 and is 
not limited to that given in Eq. 



(3, 



In the rest of 
this paper fi is set to one for typographical convenience, 
except when we discuss the semiclassical limit. 



A. Exact Solutions 

To solve the MSW equations we first introduce the 
operators H] 



A+(ai)r(ai) 



B_ =f-i(ai)i_(ai). 



(3.10) 



Using Eq. (^^) one can easily prove the commutation 
relation: 



[B^,B+]=R{ao), 
where we used the identity 

i?(a„) = r(ai)i?(a„_i)r"^(ai). 



(3.11) 



(3.12) 



valid for any n, to define Oq. One can also prove by 
induction two additional commutation relations 

[B+B^,Bl] = (i?(ai) + R{a2) + ■ ■ -fi?(a„))S!^, (3.13) 

and 

BZ"] = {R{ai) + R{a2) + • • +i?(a„))BZ" . 

(3.14) 



Using the operators introduced in Eq. (B.10),Eq. (3.4) 
can be rewritten as 



B+B-^e{x)=K^e{x). 



(3.15) 



Eqs. (3.13) and (3.14) imply that B + an d _B_ can be 



used as ladder operators to solve Eq. ( 3.15 ). To this end 



we introduce 4*^^ as the solution of the equation 
i_(ai)«'L"^ =0 = B_«'L"\ 



which implies 

*L°)^exp(-^/^(.;ar)d. 
If the function 

n 

/(n) = ^i?(afc) 
fc=i 



(3.16) 



(3.17) 



(3.18) 



can be analytically continued so that the condition 

/(Ai) = A (3.19) 



is satisfied fo r a p articular (in general, complex) va lue o f 
/i, then Eq. (3.13) implies that one solution of Eq. (3.15) 



is B'f 'I'^-'. Similarly if ^E'^"' satisfies the equation 



(3.20) 



which implies that 
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exp ( +i I ip{x; ai)dx j , 



(0) 



exp +i I ip{x; ao)dx 



(3.21) 



(3.22) 



Bl¥^^ = (^(a„+i) + ip{ai)){ipian+i) + ¥'(02)) 
••• x(.^(a„+i) + (y3(a„))exp ^-i J ip{x; an+i)dx ] , (3.28) 

and 



(0) 



(-<^(a„) - V5(ai)) ^(-(/5(a„) - ((£'(a2)) ^ 



then Eq. ( 3.14 ) implies that a second solution of Eq. 
( 3.15 ) is given by Hence for shape invariant 



electron densities the exact electron neutrino amplitude 
can be written as 



5'e(x) = PB'^exp y—i J (p(x;ai)dx 

+ 7i?Z^^^ exp I +j / ip{x;ao)dx], (3.23) 



••• x(-(/p(a„) - i^(a„)) ^ exp y+i J ip{x;an)dxj . (3.29) 

In both of these equations the quantity (^(a„) stands for 
(p{x;a„). 

If the density profile satisfies the condition 



ip{x; an) = ip{x; ai) + {n - l)C{x), 



(3.30) 



where (3 and 7 are constants to be determined by the 
initial conditions. 



B. Asymptotic Results 



then these asymptotic equations can be cast in a form 
suitable for analytic continuation : 



and 



To evaluate Eg . ( ^^2^ ) asymptotically we fir st note 
that using Eq. ( 3.10 ) the two terms in Eq. ( 3.23 ) can be 
cast in the forms 



r(n + 2(p(ai)/C) 



^-n^(0) ^ ^_^yn r(n - 1 + 2ip{ai)/() 



(3.31) 



B1 * 



(0) 



r(2n - l + 2(^(ai)/C) 
X exp [ +i / ip{x;an)dx ). (3.32) 



A+(ai)A+(a2) • •l+(a„)cxp l^-i J ip{x;an+i)dxj , (3.24) 
and 

Al^(ai)il^(a2) • •il^(a„)cxp J ip{x; an)dx^ . (3.25) 

We consider two different limits: i) f{x) is constant (i.e., 
either the electron density is constant, or the neutrino is 
traveling through the vacuum); and ii) 'p{x) —>■ ±00. In 
the former limit the commutator 



The condition in Eq. ( 3.30| ) is satisfied for a number of 
density profiles as illustrated in the next section. If this 
condition is not satisfied, the analytic continuation may 
still be done, but wil l be more complicated. Whenever 
the condition in Eq. ( 3.30| ) is satisfied, in the asymptotic 
limits the electron neutrino amplitude can be written as 



r(2^^+2^(ai)/C) , ■ I ( ^ , 

^ pC t^^.. , n.-/_ ^ /A^ ^^p -I / ip[x; a^+i)dx 



r(Ai + 2(^(ai)/C) 



^dx 



, <p(x;a„)] = ip'{x;an) 



(3.26) 



r(2/i+l + 2(^(ai)/C) 
X exp J ip{x; afi+i)dx 

Also using the identity 



vanishes. In the latter case this commutator can be ig- 
nored as 



hm l£(^.l, 



(3.33) 



(3.34) 



ip{x; a„)(/?(x; a„) + iip'{x; a„) 
= (p{x; an)(p{x; a„i) [l + i 



the asymptotic value of Eq. (3.33) for ip — > ±00 is found 
to be 



ip'{x; an) 



(p{x; an)(p{x; am) , 
(p{x;an)(p{x;arn) (3.27) 

provided that (/?'(x; a„)/(/3(x; a„) remains finite (which is 
the case for all realistic density prof iles). Hence it both 
limits we can write Eqs. ( 3.24 ) and ( 3.25| ) as 



/3(2(^(ai))''exp / (^(x; a^,+i)dx 



+ 7(-2<^(ai)) ^ "^exp^+i J ip{x;af^+i)dx^ . (3.35) 

The hopping probability can be easily obtained from 
these results as we illustrate in the next section. 
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IV. EXAMPLES 

A. Linear density profile 

For the linear density profile 

iVe(x) ^ Na- Nl^{x- xr), 
where Nq is the resonant density: 

2V2 GfNqE = dm^ cos26i«, 
we can write 

Sm'^ N' 



(4.1) 
(4.2) 

(4.3) 



Adoptin g th e ansatz 93(0;; a„) — a„(a; — xr), one finds 
that Eq.(3^) is satisfied with 



ai 



cos 29^ — = a2 = as ■ 



and 



6m? N' 
R{an) = 2iai = -i^rrr cos26'„— ^, 



(4.4) 



(4.5) 



for aU n (i.e., = for a linear density profile in Eq. 
( ^!30|) . The function /(n) of Eq. (|I|) is then 



/(") 



2E 



■ cos 29„ n 



and hence the solution of Eq.(3.19) is given by 



where we defined 



n = 



^ = 2"' 



5to2 sin2 26),„ iVo 
~AE"cos29^J^' 



(4.6) 



(4.7) 



(4.8) 



Linear density profile is special in the sense that both 
before and after the MSW resonance its absolute value 



gets to be very large. Hence Eq. (3.35) can be used to 
evaluate both asymptotic limits. Initially we take if — > 
+00 which gives 



M/,(x) =/3(2^(x))^"/' 



exp 



(p{xi)dx} 



(4.9) 



Imposing the initial condition '^{xi) = 1 determines f3 
and gives the final (tp —00) solution to be 



-I I ip{xi)dxi ]. (4.10) 



Using Eqs. ( p.5[ ) and (2.t) one finds that the asymptotic 
matter angles are 



tp — > +00, sin^i = 1, cos6'i = 0, 
if ^ —00, sinflj = 0, COS0/ = 1, 



(4.11) 



Inserting Eq. ( [l.llD into Eq. ( |2.9| ) one finds that 

^e{xf) = -n{xf). (4.12) 



Comparing Eqs. ( I.IC ) and ( 4.12 ) and using 

ip{xf) 



lim 



= -1 = 



(4.13) 



Xi^ — oo,a;/— ►+CXD ip{xi) 

we find that the hopping probability is given by 

Phop = |*2(a:/)P = exp (-ttI]) . (4.14) 

B. Exponential Density Profile 

In this case we have 

iVe(x) =iVoe""("^~^''\ (4.15) 

where Nq is the resonant density given in Eq. (^|^), we 
can write 



Adopting the ansatz 



■ cos 29. 



- 1 



(4.16) 



if{x]an) ^ ^^cos2e^{e """^ - a„) ,ai = l, 

(4.17) 

one can show that to satisfy Eq. ( |3.5D we should require 
ia'iE 



an = fli 



(Jm^ cos 29y 



(n-1), 



and 



Hence 



i?(a„) = —2iaan 



5m 



cos 26'u — 



(4.18) 



(4.19) 



(4.20) 



In this case there are two solutions of Eq. ( 3.19| ) and the 
appropriate solution giving the correct limiting behavior 
(that n vanishes as A does) is 



i ( 6m} 



a\4:E 



(1 - cos26'„). 



(4.21) 



Exponential density profile, like the linear one, can get 
to be very large (and of course positive) before the MSW 
resonance point. Hence initially the second term in Eq. 
( |3.35D vanishes and the constant p can be fixed to be /3 = 
2((^(xi))^^. However, unlike the linear density, much 
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after the MSW resonance point the exponential density 
vanishes and one should use the asymptotic form for the 
constant </? derived in the previous section. In this limit 
we have 



V(a;,a^+i) 



5m? 
"1^ 



cos 201;, 



(4.22) 



Inserting Eqs. ([4.2l[) and ( [4.22[ ) into Eq. (|3.33| ) we find 
the final electron neutrino amplitude to be 



( i Sm^ \ 
\a 2E J 



M^f)J r(i^(i + cos2(?,„)) 
1e' 



X exp +1 



+ -f{—ia) 



X exp 



r(^w(i + cos2(?„) 



r 1 



a 2B / 



,5vn? 



(4.23) 



The asymptotic matter angles are now given by 

tp +00 , sin 0i — 1, cos 6i — 0, 

95 — > 0, sin 0/ = sin , cos 6*/ = cos 01, . 

Inserting Eq. ( [4. 24) ) into Eq. ( |2.10| ) we get 

. (5m^ ^ 

L 



(4.24) 



\I/g(i) = - cos6'„\I' 

— sin 01, (-g-) exp | —i 



AE 



(4.25) 



Comparing Eqs. ( 4.23| ) and ( 4.2?: ) one gets 

1 / \ M r f i-^^ 

I / a \ ^ \a 2E ) 



2,(^) 



cos 01, \2iy9(a;/) 



X exp 



a 8E 



r 

(1 - cos 20. 



^ 4^ (l + cos20„)) 



(4.26) 



which, upon squaring, gives the hopping probability to 
be 



-Phop — 

where we defined 



g-7r(5(l-cos2e^) _ g-27r(5 



I - e 



-27r5 



6m 
2E^' 



(4.27) 



(4.28) 



Eq. ( 4.27] ) was previously obtained by direct solution 
of the differential equation The method presented 
here is not only more straightforward, but also explicitly 
demonstrates the role of the boundary conditions. Indeed 



using the infinite Landau-Zener met hod o n ly p rovides the 
first term in the numerator of Eq. ( 4.26| ) ||l5| . To obtain 
the exact result, it is essential to take into account the 
fact that the electron density goes to zero (not to infinity) 
when the neutrinos move very far away from the MSW 
resonance region as illustrated above. 



V. APPROXIMATE EXPRESSIONS 

The method presented here is useful not only to write 
down exact expressions for the electron neutrino survival 
probabilities for shape-invariant electron densities, but 
can also be used as a starting point of several useful ap- 
proximations. In particular, the transformation in Eq. 
( |3.7[ ) can be formally introduced without an explicit ref- 
erence to the parameters a„. Suppose, for example, we 
wish to satisfy Eq. ( |3.19D with n = 1 by introducing the 
function ^(z) in Eq. (|3.5|) 



^'^{x) - ifi^ix) + itvf'ix) + ih('{x) = A. 



(5.1) 



In this section we write h explicitly in the equations to 
emphasize the semiclassical nature of the approximations 
we consider. By expanding the function £^{x) in powers 
of n 



a^)^Co{x) + Hi{x) + h^U^) 



(5.2) 



we find that Eq. (5.1) can be solved by matching the 
powers of h: 



Co = + A, 



266 



(5.3) 



Using Eqs. ( |3.23| ), ( |3.24| ), and ( |3.25D the electron neu- 
trino amplitude in this approximation can be written as 

*e(a;) = /3 + exp (-^ j £,{x)dx \ + 7 



ox 



exp J ^{x)dx 



(5. 



For a monotonically decreasing cp we choose ^ 



+ y^(p^ + A. To evaluate Eq. (5.4) in the lowest order 
in h, we note that in the exponentials one needs to in- 
clude ^1 (since ?i's cancel) whereas in the pre-exponential 
factors it is sufficient to retain only ^q. Inserting Eq. 
( ^.3[ ) into Eq. (5.4) and imposing the initial conditions 
we obtain to lowest order in h 

"feix) = ir_(0)T_(a:)exp (^+^ J ^p^x)+Adx^ 

+ ir+(0)r+(x)exp('-J / ^p^x)+Adx) , (5.5) 



where we defined 
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1 ± 



1/2 



Eq. (5.5) is the electron neutrino amplitude in the stan- 
dard adiabatic approximation. The equivalence of the 
adiabatic and primitive semiclassical approximations was 
previously proved in Ref. Here we demonstrated 

that it can be obtained by solving a shape-invariance con- 
dition, valid for any electron density, to lowest order in 

n. 

Even though shape-invariant electron densities consti- 
tute a rather restrictive class, one can use these exact so- 
lutions as comparison solutions to obtain uniform semi- 
classical approximate solutions [ p^UTl|Jl7[| for any given 
density that is not shape- invariant. To achieve this goal, 
one starts with a shape-invariant electron density profile 
for which the solution of the "mapping" equation 



matter-enhanced neutrino oscillations, namely linear and 
(5.6) exponential densities. In addition to these two the tanh a; 

and l/{a + hx) density profiles are also shape-invariant 
and the MSW equations for these density profiles can be 
solved using the techniques described here. One should 
emphasize that although the neutrino propagation in 
these density profiles was studied before, the method de- 
scribed in this paper seems to be more direct and to have 
pedagogical value as it explicitly illustrates the role of the 
boundary conditions. The resulting neutrino amplitudes 
can also be utilized as comparison amplitudes for the uni- 
form semiclassical treatment of neutrino propagation in 
arbitrary electron density profiles. 



d_ 

dS 



as) xiS) - nxiS) (5.7) 



5_2) 



is known. In the next step, the variable S in Eq 
is taken to a function of x and the solutions of Eq. (3.4) 
are written as 



<f{x)=K{x)xiSix)). 
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K{x) = 



(5.9) 



one can show that ^ given in Eq. (5.8) satisfies Eq. (3.4) 
if the equality 



K" 



K 



+ (A + (^^ + ihip') = 



(5.10) 



is satisfied. Eq. ( ^.10| ) is usually solved by expanding 
S{x) in powers of h. Using a linear mapping density 
this mapping technique was used in Refs. |1^ and p^ ] 
to obtain an approximate solution for the neutrino sur- 
vival probability. It is possible to generalize this result by 
using not the linear density, but another shape-invariant 
mapping density, which most closely resembles the actual 
electron density that the neutrinos travel through. 



VI. CONCLUSIONS 

We showed that the analogy between supersymmetric 
quantum mechanics and matter-enhanced neutrino os- 
cillations can be used to obtain exact expressions for 
the neutrino amplitudes traveling in a class of electron 
density profiles. We explicitly worked out the neutrino 
survival probability for two shape- invariant electron den- 
sity profiles that are most commonly used for studying 
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